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We report on the results of penguin-mediated B decays at the Belle experiment. The analyses 
were based on approximately 32 million BB events collected at the T(4S) resonance with 
the Belle detector at the KEKB e + e~ storage ring. The b — ► sj transition was studied 
through exclusive decays: B -> K*^, B° -> A" 2 *(1430)°7, B + -> K+^-K'-y, B + -> K*°n + 'y 
and B+ -> K+p°f. The b -> si + £~ transition was searched through both exclusive decays, 
B — > K^t + l~ , and inclusive decay, B — > We observed the decay processes £? + — > 

_ft' + 7r + 7r~7 and _B —* K£ + £~ for the first time. 



1 Introduction 



In the Standard Model (SM), flavor-changing neutral current (FCNC) decays are forbidden 
at tree level. However, FCNC decays are induced through loop diagrams, such as penguin 
diagrams or box diagrams. These loop diagrams are sensitive to new physics, since heavy 
particles beyond the SM, such as charged Higgs or SUSY particles, contribute to additional loop 
diagrams, branching fraction or kinematic variables can be deviated from the SM values. 

We have studied the FCNC decays of b — > S7 and b — ► sl + i~ using a data sample collected 
with the Belle detector 111 at the KEKB El storage ring. The data sample corresponds to 29 fb 
taken at the T(45) resonance and contains approximately 32 million BB pairs. 



2 b — > 57 transition 

2.1 Analysis of B — > K*^ 

Exclusive B — > K*^ decays were reconstructed from a high energy photon and a K* (K*° — > 
K + it~ or Ks^°, K* + — > Ksn + or K + tt°). The continuum background was suppressed by a 
likelihood ratio, which was constructed from an event shape variable Super Fox- Wolfram (SFW) 



and the B meson flight direction. After applying a likelihood ratio cut, the K*"f final state 
was cleanly reconstructed. The beam energy constrained mass (Mb c ) distributions for each sub 
decay modes are shown in Fig p]. We obtained the branching fractions: 

B(B° K*°j) = (4.08j;°;! ± 0.26) x 10~ 5 , (1) 
B(B+^K*+ 7 ) = (4.92±g;||±°;if) x HT 5 . (2) 

We also checked for any partial decay rate asymmetry. We only used the self-tagging modes: 
-fT ± 7r^7, Er ± 7r°7 and Kstt ^y. The wrong tag fraction due to hadron misidentification was 
estimated to be only 1.2%. We determined the partial rate asymmetry as 

A CP = n-B^-TiB^K*,) = 

F(B — > K 7) + T(B — > K*~f) V ~ 6 " 8 ; KJ 

which corresponds to 

(-8.5 < A C p < 14.9)% (90% C.L.), (4) 

and was found to be consistent with zero. This limit is the most stringent over previously 
published results §§,§]]. 



40 



20 



g 10 



PRELIMINARY 



(a) KVy 



Yield 180.5 +15.5-14.8 +2.7-3.1 




(C) K + 7t°Y 




(b) K S 7T + Y 



Yield 64.2 +9.4-8 7 +1.2-1.1 




Yield 7.3 +3.4-2.8 ± 0.2 

fit result 

total b.g. 



>.2 5.22 5.24 5.26 5.28 5.2 5.22 5.24 5.26 5.28 

Beam constrained mass (GeV/c 2 ) 



20 
15 
10 
5 

> 

°t 

4 .9 



- 3 



Figure 1: M bc distribution of B -> K* (892)7- 



j8.j8 Analysis of B° -> K|(1430)°7 

In the -B — > iT|(1430)°7 analysis, the selection criteria used were similar to those used in 
the if* 7 analysis. The if|(1430)° was reconstructed from K + ir~, which lies between 1.25 and 
1.60 GeV/c 2 . We obtained 27.0 ± 6.7±§'| events for the B -» if+TT^ decay from the M bc fit 
(Fig. 1). 

To distinguish the signal B° -> if|(1430)°7 from Ki (1430)7 and non-resonant iT + 7T 7, we 
performed a unbinned maximum likelihood fit to Mb c ,the helicity angle (cos#hei) and the K + n~ 



invariant mass (Mk-k)- Fig. Jl shows the cos#hei distribution, where the continuum background 
has been subtracted. We found that -fT|(1430) resonance was the dominant component and 
obtained a signal yield of 24.0^§ 4 events. The branching fraction of B° — > iT|(1430) 7 was 
determined to be 

B(B° - K 2 *(1430)° 7 ) = (1.50±8:g±8:i|) x 10- 5 . (5) 
This result is consistent with a previous measurement H and theoritical predictions.! 
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Figure 2: M^ c distribution (left) and cos#hci distribution (right) in B — > K + 7v 7 decay. 



2.3 Analysis of B + — ► K + tt + tt^'j 

We extended the same analysis to the three body hadronic final states of radiative decay. The B + 
was reconstructed from K + ir + ir~ r y. The invariant mass of K + ir + n~ was required to be between 
1.0 GeV/c 2 and 2.4 GeV/c 2 . To extract the signal yield, we made a fit to the Mb c distribution 
and obtained ST.TtJx'ili'g signal events with a statistical significance of 6. Oct (Fig. |^). We first 
observed B + — ► ET + 7r + 7r - 7 and determined the branching fraction as 

B(B+ - K+tt+tt-t) = (2.43t°;f 7 t°;l) * (6) 

The invariant mass of reconstructed K +, k + it~ is shown in Fig. ||. We can observe a significant 




Figure 3: The Mb c distribution of B — > Kirny. 

excess below 1.8 GeV/c 2 . In this mass region, there are many higher kaonic resonances which 
contribute to this decay mode, which are difficult to distinguish. However we can still identify 
K + ir + ir~ via K*ir or Kp. We performed on unbinned maximum likelihood fit to M^ c , the 



invariant mass of K + tt~ (Mktt) and 7r + 7r - (M^). Fig. |] shows the Mk-k and M n7T distributions 
along with the fit. We obtained the signal yield of K*tt^, Kp-y and non-resonant Kirirj to be 
32.5^j;o'o-i'7) 24.2 ± H.6I7Q and O.OlJ 1 ^ , respectively. The branching fraction and upper limits 
were found to be 

B{B + -» iTOvr+7) = 2.04l°;^ ± 0.22) x 10~ 5 , (7) 
B(B + -> K+p°j) < 1.9 x 10~ 5 (90% C.L.), (8) 
B(B + -> i^ + vr + vr"7 non - resonant) < 0.92 x 10" 5 (90% C.L.). (9) 

This result shows that the B + — > lf + 7r + 7r~7 process is consistent with a mixture of B + — > 
K*0vr+7 and 5+ -> K+p ^. 

We measured the exclusive two and three body hadronic final states of radiative decay. 
These exclusive decay rates could be compared with the inclusive B — > X s 7 decay rate.H When 
we calcurate the exclusive decay rate, the isospin invariance in the decay rate was assumed. 
We summed up those decay rates, giving 11.7 ± 2.0 x 10~ 5 . This accounts for 35 ± 11% of the 
inclusive decay rate. 




Figure 4: The Mk-k (left) and M nn (right) distributions in B ~ * Knnj decay. 



3 b — * s£ + £ transition 

3.1 Analysis of B —> K^£+£- 

Here, we summarize what is already published in Ref.i We reconstructed B — > K^£ + £~ from 
oppositely charged lepton pair (e + e _ or p + p~) and a kaon or a K* (K + , K$ — > 7r + 7r~, K*° — > 
K + 7T~ or ^ 5 7r , -» K 5 vr+ or K+7T ). The JytyX, and ^ X s events were vetoed by di- 

lepton invariant mass cuts. The di-electron which comes from photon conversion and 7r° Dalitz 
decay are removed by requiring M ee > 0.14 GeV/c 2 . The continuum background was suppressed 
using a likelihood ratio formed by a Fisher discriminant, cos 9b and the angle between the B 
candidate sphericity axis and z axis cos # S ph- The Fisher discriminant was calculated from the 
energy flow in 9 cones along the B candidate sphericity axis and the normalized second Fox- 
Wolfram moment (-ffe)- Another major background is BB events where both B mesons decay 
to X c lv. This background is suppressed by another likelihood ratio constructed from a missing 
energy in the event and cos 8b- 

To extract the signal yield, we made a fit to the Mb c distributions (Fig. ||). We observed 
Q.Sltl'ilti events for B — > Kp^ pj~ with a statistical significance of 4.7. The di-muon invariant 
mass distribution is shown in Fig. ^, and is found to be consistent with the Monte Carlo (MC) 
expectation. For other modes, we observed no significant excess: 4.1^2 1^08 events for B — > 



Ke + e~, 6.3±|;5±i;? events for B -» i^e+e" and 2.1±|;f±?5 events for B -» K*n + fjT. In the 
combined K£ + £~ mode, we observed 13.6l 3 [ 8 lij events with a statistical significance of 5.3. 
For the mode with significant signal events, we determined the branching ratios: 

B(B^K£ + r) = (0.7518;!? ±0.09) x 1(T 6 , (10) 
B(B — > K[i + fi~) = (0.99l°;f°l°il) x 10" 6 , (11) 

and found them to be consistent with the SM predictions.il For other modes, we set the upper 
limits as follows: 

B(B -> Ke + e~) < 1.3 x 10~ 6 (90% C.L.), (12) 
B(B -» iTe+e") < 5.6 x 10~ 6 (90% C.L.), (13) 
B(B -> K*ii + fi-) < 3.1 x 10" 6 (90% C.L.). (14) 
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Figure 5: Afbc distributions of B — » 7^'*^ + ^ _ decay. The left column is for the data and the right column is for 

the MC background. 
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Figure 6: M w distribution of B — > K*£ + £ . The hatched histogram shows the data distribution, while the open 

histogram shows the MC signal distribution. 



3.2 Analysis of B —> X s i + £~ 



The B — > X s ^ + -£ _ decay was reconstructed by combining an X s system and an oppositely 
charged lepton pair ( e + e~ or /i + ). The X s is formed from one charged or neutral kaon 
and zero to four pions where at most one neutral pion is allowed. The continuum background 
is suppressed by R2 and the he remaining tracks. The BB background is suppressed by two 
likelihood ratios. The first likelihood ratio is formed from the missing energy and the X s mass. 
The second likelihood is constructed from cos 6b and the sum of the cosine of the angle between 
the kaon and leptons {cos9 K i+ + cos6 K £-). The best candidate is selected by AE. 

The signal yield is extracted from a fit to the Mb c distributions. We have found no significant 
excess, and set the upper limits as follows: 

B(B^X s e + e~) < 10.2 x 1(T 6 (90% C.L.), (15) 
B{B -> Xsn + n~) < 19.9 x 10~ 6 (90% C.L.). (16) 

These results are close to the SM predictions, i 



4 Summary 

We have studied penguin-mediated B decays. The branching fractions of B — * K*~f, B° — > 
-K"!(1430)°7, B + — s> iT + 7r + 7r~7, B + — > K*°tt + 'j were measured and the upper limit of B° —> 
7 was set. The B + — » K + tt + tt 7 was observed for the first time. In B — ► K*j decay, we 
set the most stringent limit on the partial rate asymmetry, which was found to be consistent 
with zero. 

We observed B — > K£ + £~ , and the measured branching fraction was found to be consistent 
with the SM predictions, i For other b — » s£ + £~ transitions, we set the upper limits of the 
branching fractions. These results were used to constrain the Wilson coefficients Cg ff and Ciq." 
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